Influence of strain on anisotropic thermoelectric transport of Bi2Te3 and Sb2Te3 
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On the basis of detailed first-principles calculations and semi-classical Boltzmann transport, the 
anisotropic thermoelectric transport properties of Bi2Te3 and Sb2Te3 under strain were investigated. 
It was found that due to compensation effects of the strain dependent thermopower and electrical 
conductivity, the related powerfactor will decrease under applied in-plane strain for Bi2Te3, while 
being stable for Sb2Te3. A clear preference for thermoelectric transport under hole-doping, as well 
as for the in-plane transport direction was found for both tellurides. In contrast to the electrical 
conductivity anisotropy, the anisotropy of the thermopower was almost robust under applied strain. 
The assumption of an anisotropic relaxation time for Bi2Te3 suggests, that already in the single 
crystalline system strong anisotropic scattering effects should play a role. 

PACS numbers: 31.15.A-,71.15.Mb,72.20.Pa,72.20.-i 



I. INTRODUCTION 

Thermoelectric (TE) materials are used as solid state 
energy devices which convert waste heat into electric- 
ity or electrical power directly into cooling or heating^^'^. 
TcUuridc based thermoclcctrics, e.g. the bulk materials 
bismuth (Bi2Te3) and antimony telluride (Sb2Te3) and 
their related alloys, dominate efficient TE energy conver- 
sion at room temperature for the last 60 years'''^. The 
materials TE efficiency is quantified by the figure of merit 

ZT = T, (1) 

where a is the electrical conductivity, S the thermopower, 
Kei and are the electronic and phononic contribution 
to the thermal conductivity, respectively. From Eq. 1 
it is obvious, that a higher ZT is obtained by decreas- 
ing the denominator or by increasing the numerator, the 
latter being called powerfactor PF — crS^. While bulk 
Bi2Te3 and Sb2Te3 show values smaller 1 and appli- 
cations have been limited to niche areas, a break-trough 
experiment of Vcnkatasubramanian et al. showed a re- 
markable ZT = 2.4/1.5 for p-typc/n-type superlatticcs 
(SL) composed of the two bulk tcUurides^"^. With the 
availability of high-ZT materials, many new applications 
will emerge'^. The idea of thermoelectric SL follows the 
idea of phonon-blocking and electron-transmitting at the 
same time. It suggests that cross-plane transport along 
the direction perpendicular to the artificial interfaces of 
the SL reduces phonon heat conduction while maintain- 
ing or even enhancing the electron transport"^. While 
some effort in experimental research was done^"-'^'^, only 
a few theoretical works discuss the possible transport 
across such SL structures^^'^^. While Park et al}'^ dis- 
cussed the effect of volume change on the in-plane ther- 
moelectric transport properties of Bi2Te3, Sb2Te3 and 
their related compound, Li et al}^ focussed on the cal- 
culation of the electronic structure for a Bi2Te3/Sb2Te3- 
SL, stating changes of the mobility anisotropy estimated 
from effective masses. 



Superlattices are anisotropic by definition and even the 
telluride bulk materials show intrinsic anisotropic struc- 
tural and electronic properties. However, investigations 
of Vcnkatasubramanian et al. found a strong decrease 
for the mobility anisotropy and the thermoelectric prop- 
erties for the Bi2Te3/Sb2Te3-SL at certain periods. The 
reason for this behaviour is still on debate and could be 
related to strain effects which are induced by the epitax- 
ial growth of the Bi2Tc3/Sb2Te3-SL. To extend previous 
works^^"^^ and to clarify the open question on the re- 
duced anisotropy, we are going to discuss in this paper 
the anisotropic electronic transport in bulk Bi2Te3 and 
Sb2Te3 and the possible influence of strain in epitaxially 
grown SL on the TE properties. 

For this purpose the paper will be organized as follows. 
In section II we introduce our first principle electronic 
structure calculations based on density functional the- 
ory and the semi-classical transport calculations based on 
the solution of the linearized Boltzmann equation. With 
this, we discuss the thermoelectric transport properties, 
that is electrical conductivity, thermopower and the re- 
lated powerfactor, of unstrained Bi2Te3 and Sb2Te3 with 
a focus on their directional anisotropics. While in epi- 
taxially grown Bi2Te3/Sb2Te3-SL the atoms near the in- 
terfaces may be shifted from their bulk positions due 
to the lattice mismatch and the changed local environ- 
ment, we modelled Bi2Te3 with the experimental lat- 
tice parameters and interatomic distances of Sb2Te3 , 
and vice versa. We assume that from these two limit- 
ing cases one could estimate the effect of the interface 
relaxation on the electronic and transport properties in 
Bi2Te3/Sb2Te3-SL. With that structural data we first 
analyse in section III the anisotropic thermoelectric prop- 
erties of the unstrained bulk systems, while in section IV 
a detailed view on the influence of strain, which may oc- 
cur in Bi2Te3/Sb2Te3-SL, on the electronic transport of 
these tellurides is given. Throughout the paper we quote 
Bi2Te3 (Sb2Te3) as strained, if it is considered in the lat- 
tice structure of Sb2Te3 (Bi2Te3). As in the SL p-type. 
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as well as n-type, transport was reported, we discuss the 
concentration dependence for both types of carriers on 
the transport properties. 



II. 



METHODOLOGY 



For both bismuth and antimony telluridc we used 
the experimental lattice parameters and relaxed atomic 
positions^^ as provided for the rhombohedral crystal 
structure with five atoms, i.e. one formula unit, per unit 
cell belonging to the space group Z?|^ (i?3m) . The related 
layered hexagonal structure is composed out of three 
formula units and has the lattice parameters a!^^^ = 
4.384A, c^-J^g = 30.487A, and a^|f,^ = 4.264A, cf^^ = 
30.458A, for Bi2Te3 and Sb2Te3, respectively. In fact, 
the main difference between the lattices of Bi2Te3 and 
Sb2Te3 is a decrease of the in-plane lattice constant with 
an accompanied decrease in cell volume. So, a change 
between the two lattice constants can be related to ei- 
ther compressive or tensile in-plane strain. This is very 
similiar to the approach by Park et al}^ , while omitting 
computational relaxation of internal atomic positions. 

Our electronic structure calculations arc performed in 
two steps. In a first step the detailed band structure 
of the strained and unstrained Bi2Te3 ans Sb2Te3 was 
obtained by first principles density functional theory cal- 
culations (DFT), as implemented in the fully relativis- 
tic screened Korringa-Kohn-Rostoker Greens-function 
method (KKR)^°. Within this approach the DiRAC- 
equation is solved self-consistently and with that spin- 
orbit-coupling is included. Exchange and correlation ef- 
fects were accounted for by the local density approxima- 
tion (LDA) parametrized by Vosco, Wilk, and Nusair^^. 
A detailed discussion on the influence of strain on the 
band structure topology of Bi2Te3 and Sb2Te3 is recently 
published^^. With the well converged results from the 
first step we obtain the thermoelectric transport prop- 
erties by solving the linearized Boltzmann equation in 
relaxation time approximation (RTA) within an in-house 
developed Boltzmann transport code^'^"^^. Boltzmann 
transport calculations for thermoelectrics have been car- 
ried out for quite a long time and show reliable results 
for metals^^"^^ as well as for wide- and narrow gap semi- 
conductors^^'^^"'^^. TE transport calculations for bulk 
Bi2Te3"A8,33,34 Sb2Te3"^i^'35 ^gj.g presented be- 
fore. Here the relaxation time t is assumed to be con- 
stant with respect to wave vector k and energy on the 
scale of fc^T. This assumption is widely accepted for 
metals and highly doped semiconductors. Most of the 
presented results are in this high-doping regime. Within 
the RTA the transport distribution function (/i, 0) 
(TDF) and with this the generalized conductance mo- 
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FIG. 1: (color online) Band structures of (a) Bi2Te3 and (b) 
Sb2Te3 along symmetry lines for both unstrained (black solid 
lines) and strained (red dashed lines) lattices. Energies are 
given relative to the valence band maximum. 

ments -C^ji (/i, T) are defined as 



E=EY. 



(2) 

^fc (II)' ''^k (_L) denote the group velocities in the directions 
in the hexagonal basal plane and perpendicular to it. 
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respectively. Within here the group velocities were ob- 
tained as derivatives along the lines of the Blochl mesh 
in the whole Brillouin zone^^. A detailed discussion on 
implications and difficulties on the numerical determina- 
tion of the group velocities in highly anisotropic mate- 
rials is currently published elsewhere'^^. As can be seen 
straight forwardly, the electrical conductivity a in- and 
cross-plane is then given by 



(3) 



and the temperature- and doping-depcndcnt ther- 
mopower states as 



1 4'VM,r) 



(4) 



for given chemical potential /i at temperature T and ex- 
trinsic carrier concentration N determined by an integra- 
tion over the density of states niE) 

where CB""" is the conduction band minimum and 
ygmax valence band maximum. The energy range 

AE has to be taken sufficiently large to cover the tails 
of the Fermi-Dirac distribution function f(fj.^T) and to 
ensure convergence of the integrals in eq. 2 and 5^^ . The 
k-space integration of eq. 2 for a system with an intrin- 
sic anisotropic texture is quite demanding. In previous 
publications^^''^'' we stated on the relevance of adaptive 
integration methods needed to reach convergence of the 
energy dependent TDF. Especially in regions close to 
the band edges the anisotropy of the TDF requires a 
high density of the k-mesh. Here, convergence tests for 
the transport properties showed that at least 150 000 k- 
points in the entire BZ had to be included for sufficient 
high doping rates {N > 1 x 10^^ cm~^), while for energies 
near the band edges even more than 56 million k-points 
were required to reach the analytical values for the con- 
ductivity anisotropies at the band edges^'^. Within the 
RTA, from comparison of the calculated electrical con- 
ductivities (eq. 3) with experiment it is possible to con- 
clude on the directional anisotropy of r. For the ther- 
mopowcr S (eq. 4) the dependence of the TDF on the en- 
ergy is essential. That is, not only the sloop of the TDF, 
moreover the overall functional behaviour of the TDF 
on the considered energy scale has to change to observe 
an impact on the thermopower. The calculations in this 
paper aim to cover band structure effects and not scatter- 
ing specific impacts by an energy- and state-dependent 
relaxation time. 
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FIG. 2: (color online) Anisotropic thermopower for bulk (a) 
Bi2Te3 and (b) Sb2Te3 in their unstrained bulk lattice con- 
stants. Electron doping refers to the blue (thick) lines in 
the lower part of the figure, while red (thin) lines refer to 
hole doping and positive values of the thermopower. Solid 
lines show the in-plane part S^^ of the thermopower, while 
dashed lines show the cross-plane part S ^ . The extrin- 
sic charge carrier concentration of Bi2Te3 and Sb2Te3 was 
fixed to A'' = 1 X 10" cm"^ and = 1 x 10^° cm"^, respec- 
tively. Experimental data (squares, diamonds, circles, trian- 
gles) from Ref. 41-43 are given for comparison. 



III. ANISOTROPIC THERMOELECTRIC 
PROPERTIES OF UNSTRAINED Bi2Te3 AND 
Sb2Te3 

In order to understand the experimental find- 
ings on the in-plane and cross-plane transport of 
the Bi2Te3/Sb2Te3-SL, in the following section the 
anisotropies of the electrical conductivity, the ther- 
mopower and the related powerfactor of bulk Bi2Te3 and 
Sb2Te3 are discussed. Even though the behaviour of 
Sb2Te3 is strongly p-typc with a extrinsic carrier con- 
centration of iV = 1 . . . 10 x lO'^^ cm^^,-^^ we also discuss 
the related n-doped case, as in Bi2Tc3/Sb2Tc3-SL n- as 
well as p-doping was reported. Bulk Bi2Te3 is known to 
be inherent electron conducting, while hole doping is ex- 
perimentally achievable for bulk systems^''^^^'*^. Figure 
2 shows the variation of the anisotropic thermopower 
for unstrained Bi2Te3 and Sb2Te3 in a wide tempera- 
ture range. The extrinsic charge carrier concentration 
of Bi2Te3 and Sb2Tc3 was fixed to iV 1 x lO^^ cm-^ 
and = 1 X 10^° cm""^, respectively. As a reference 
experimental values for both single crystalline materials 
at the same doping conditions are given and an excel- 
lent agreement can be stated. It is worth noting, that 
within eq. 4 the calculation of the thermopower is com- 
pletely free of parameters. For Bi2Te3 the in-plane ther- 
mopower reaches a maximum of S^^ — 200/xV/K at 
300 K, while the maximum for the hole-doped case is 
shifted to slightly higher temperatures of 350 K with a 
maximum values of S*,, ~ 225 //V/K. We note, that the 
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temperature of the maximum is slightly overestimated. 
This might be caused by the missing temperature de- 
pendence of the energy gap, which was determined as 
Eg = 105 meV for unstrained Bi2Te3. The anisotropy 
of the thermopower is more pronounced for the p-doped 
case. Here the cross-plane thermopower S*^ is for the 
given doping always larger than the in-planc part 6*11 . 
The anisotropy ^n/s^ is about 0.64 at lOOK, evolving to 
SJs^ - 0.79 and SJs^ - 0.55 at 300K and 500K, re- 
spectively. The sole available experimental data show no 
noticeable anisotropy for the thermopower in the hole- 
doped case'*^. For the electron-doped case the situa- 
tion is more sophisticated. While upto 340K the overall 
anisotropy is rather small, with values ^\\/s^ 0.9, a 
considerable decrease of Sj^ at higher temperatures leads 
to high values of ^n/s^ for temperatures above 400K. 
This tendency could also be revealed by experiments^^'^^. 
The crossing point of S^^ and Sj_ near room temperature 
could explain the fact of varying measured anisotropies 
for the thermopower at 300K. Here anisotropy ratios of 
= 0.97... 1.10 were reported^-^'^^. The maximum 
peak of the thermopower near room temperature can be 
explained by the position of the chemical potential as a 
function of temperature at a fixed carrier concentration. 
For T much smaller than 300K the chemical potential is 
located in either the conduction- or valence band with 
the tails of the Fermi-Dirac-distribution in eq. 2 only 
playing a subsidiary role. For rising temperatures the 
chemical potential shifts towards the band edges and S 
maximizes. At these conditions the conduction is mainly 
unipolar. For higher temperatures the chemical potential 
shifts into the bandgap and conduction becomes bipo- 
lar leading to a reduced thermopower. For the case of 
Sb2Te3 , shown in fig. 2(b), the situation is different. 
Due to the ten times higher inherent doping and the 
smaller energy gap of Eg = 90 meV, the chemical po- 
tential is located deeply in the bands for the whole rel- 
evant temperature range. Therefore the functional be- 
haviour can be understood in terms of the well known 
MOTT relation, where equation 4 qualitatively coincides 
with S (xT ■ '^''j'^^"^'' \E=fi for the thermopower in RTA^^. 
With increasing temperature the thermopower increases 
almost linearly, showing values of S^^ ^ 87/^V/K and 

~ — 72/iV/K at 300K for p- and n-doping, respec- 
tively. The anisotropy of the thermopower for the hole- 
doped case is around ^\\/s^ = 0.91, almost temperature- 
independent and slightly underestimates the available ex- 
perimental values'^ While for the electron-doped case 
the absolute values of the in-plane thermopower are com- 
parable to those of the hole-doped case, the anisotropies 
are rather large. The anisotropy varies only weakly on 
temperature showing ^\\/Sj_ = 0.48 . . . 0.52 over the hole 
temperature range. While bulk Sb2Te3 states a strong 
p-character due to inherent defects, we note here again, 
that n-doping is available in hcterostructures combining 
Bi2Te3 and Sb2Te3^ 

A strongly enhanced cross-plane thermopower 

could lead to a strongly enhanced powerfactor PF^ , 
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FIG. 3: (color online) Ratio a^^ /a^ of the electrical con- 
ductivites at 300K for unstrained bulk (a) Bi2Te3 and (b) 
Sb2Te3 . Electron doping refers to blue lines, while red 
lines refer to hole doping. The dashed lines in panel (a) 
present the ratio obtained with an anisotropic relaxation time 
txx/t-j — 0.47, while all other results are obtained with an 
isotropic relaxation time. Experimental data (circles and tri- 
angles) from Ref. 39,48,49 are given for comparison. 



if the cross-plane electrical conductivity is maintained 
at the bulk value. For this purpose the anisotropy of the 
electrical conductivity in dependence on the in-plane con- 
ductivity (T|| for unstrained Bi2Te3 and Sb2Te3 is shown 
in Figure 3. The temperature is fixed at 300K, blue and 
red lines refer to electron- and hole-doping, respectively. 
From comparison with experimental data^^, the in-plane 
relaxation time is determined to be ry = 1.1 x 10~^^ s and 
T|[ = 1.2 X 10~^^s for Bi2Te3 and Sb2Te3 , respectively. 
With that we find strong anisotropics for the electrical 
conductivity "'i\/a± ^ 1, clearly preferring the in-plane 
transport in both bulk tellurides. For the strongly 
suppressed cross-plane conduction p-type conduction is 
more favoured than n-type conduction. For Bi2Te3 the 
pure band structure effects (solid lines in Figure 3(a)) 
overestimate the measured anisotropy ratio'^^ of the 
electrical conductivity. With an assumed anisotropy of 
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the relaxation time of ■^ii/t^ = 0.47 the experimental 
values are reproduced very well. That means, scatter- 
ing effects strongly affect the transport and electrons 
travelling along the basal plane direction are scattered 
stronger than electrons travelling perpendicular to the 
basal plane. The origin of this assumed anisotropy has 
to be examined by defect calculations and resulting 
microscopic transition probabilities and state dependent 
mean free path vectors. It is well known, that in 
Bi2Te3 mainly anti-site defects lead to the inherent con- 
duction behaviour"^^'^^'^^. We have shown elsewhere'^'', 
that the integration of the transport integrals eq. 2 in 
anisotropic k-space requires large numeric effort. Tiny 
regions in k-space close to the band gap have to be 
scanned very carefully and the texture in k-space has a 
drastic influence on the obtained anisotropy values, if 
integrals are not converged with respect to the k-point 
density. As shown, some integration methods tend for 
the given k-space symmetry to underestimate the ratio 
o'li/cr^ in a systematic manner and therefore would shift 
anisotropy closer to the experimental observed values, 
without representing the real band structure effects. For 
unstrained Bi2Te3 the electrical conductivity anisotropy 
is highest for low values of a^^ , i.e. small amounts 
of doping and bipolar conduction. For larger charge 
carrier concentrations, i.e. the chemical potential shifts 
deeper into either conduction or valence band, the 
in-plane conductivity a^^ increases and the ratio '^ii/o-^ 
decreases. Values for ■^11/0-^ will lower from 7 to 2 for 
p-type conduction and 9 to 3 for n-type conduction. 
However, cross-plane electrical transport is always more 
suppressed for n-type carrier conduction, which also 
holds for unstrained Sb2Te3 . As shown in Figure 3(b) 
'^ii/tTi is almost doping independent for hole-doping, 
showing an anisotropy of around 2.7 in very good 
agreement with experiment (circle and triangles in fig. 3 
from Ref. 47-49). In this case no anisotropic relaxation 
times had to be assumed. For electron doping the ratio 
'^ii/a-± is clearly higher, evolving values of 3.5 to 6 for 
rising in-plane conductivity. The dependence of the 
anisotropy ratio on the applied doping, i.e. changing 
CTii , can be directly linked to the functional behaviour of 
the TDF near band edges, which is crucially influenced 
by the topology of the band structure^^. 



IV. ANISOTROPIC THERMOELECTRIC 
PROPERTIES OF STRAINED Bi2Te3 AND SbgTes 

Before the influence of in-plane strain on the resulting 
powerfactor will be discussed, we will first note on the 
strain induced changes of the components electrical con- 
ductivity and thermopower. In Figure 4 the anisotropy 
of the electrical conductivity '^11/0-^ is shown for both 
Bi2Te3 in the lattice constant of Sb2Te3 , i.e. under 
biaxial compressive in-plane strain (Figure 4(a)), and 
Sb2Te3 in the lattice constant of Bi2Te3, i.e. under biax- 
ial tensile in-plane strain (Figure 4(b)). 
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FIG. 4: (color online) Conductivity ratio a^^ ja ^ of the 
electrical conductivites at 300K for bulk (a) Bi2Te3 in the 
Sb2Te3 structure and (b) Sb2Te3 in the Bi2Te3 structure. 
Electron doping refers to blue lines, while red lines refer to 



and r = 1.2 x 10 s for and 
Bi2Te3 and Sb2Te3 , respectively. 



are assumed for 



For Bi2Te3 the compressive in-planc strain causes an 
increase of the the band gap by around 23% yielding 
Eg ~ 129 meV. While the anisotropy '^11/0-^ for hole dop- 
ing (red lines in fig. 4(a)) decreases to around 4 and 
is almost constant under varying doping level, the ra- 
tio considerably raises under electron doping to values 
up to 13 for tj^i ~ 100 . . . 1000 (r2 cm) — 1, correspond- 
ing to electron charge carrier concentrations of iV = 
3 ... 30 X 10^^ cm~'^. This concludes, that the cross-plane 
electrical conductivity of Bi2Te3 under compressive in- 
plane strain will be noticeably enhanced for p-doping, 
but drastically suppressed for n-doping. Such a com- 
pressive in-plane strain could be introduced by either 
a substrate with smaller in-plane lattice constant, e.g. 
GaAs-[lll] with a = 3.997A, or a considerable amount 
of Sb2Te3 in the Bi2Te3/Sb2Te3-SL. For tensile in-plane 
strained Sb2Te3 the impact on the electrical conductivity 
ratio '^ii/o-^is less prominent. As shown in figure 4(b) at 
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FIG. 5: (color online) In-plane (solid lines) and cross- 
plane (dashed lines) doping-dependent thermopower at 300K 
for (a) Bi2Te3 in the Bi2Te3 structure, (b) Bi2Te3 in the 
Sb2Te3 structure, (c) Sb2Te3 in the Bi2Te3 structure and (d) 
Sb2Te3 in the Sb2Te3 structure. Electron (hole) doping is 
presented as blue thick (red thin) line. The black (dashed- 
dotted) line in panel (d) shows the PiSARENKO-dependence of 
the thermopower expected for parabolic bands. Experimen- 
tal data (circles) from Ref. 51 is given for comparison. The 
charge carrier concentration is stated in units of e/uc (1/cm'') 
at the bottom (top) x-axis. 



hole doping '^ii/(tj_ ~ 2.5 is only marginally altered com- 
pared to the unstrained case (comp. fig. 3(b)). Mean- 
while '^ii/(7i decreases noticeably for n-type doping yield- 
ing about 3 at low a^^ and low electron charge carrier 
concentrations, and slightly higher values of '^ii/a± ~ 4 
for higher doping. Overall, the tensile strain reduces the 
electrical conductivity anisotropy by a factor of about 
1.5, directly leading to an enhanced electrical conduc- 
tivity along the z-axis of single crystal Sb2Tc3 . Wc 
note, that tensile in-plane strain opens the gap remark- 
ably by around 56% compared to the unstrained case 
to a value of Eg = 140 meV. Furthermore, such ten- 
sile strain could be incorporated by using either a sub- 
strate with larger in-plane lattice constant, e.g. PbTe- 
[111] with a = 4.567A, or a higher fractional amount of 
BiaTca in BiaTcg/SbaTea-SL. In Figure 5(a), (d) ((b). 
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FIG. 6: (color online) Change of the in-plane thermopower 
^ll under applied strain for (a) Bi2Te3 and (b) Sb2Te3 . Given 
is the ratio of in the "smaller" lattice of Sb2Te3 divided 
by in the "larger" lattice of Bi2Te3 . The doping was fixed 
to iV = 1 X 10^^ cm"^ for Bi2Te3 and iV = 1 x 10^° cm"^ for 
Sb2Te3 . Solid blue (dashed red) lines refer to electron (hole) 
doping, respectively. 



(c)) the doping dependent anisotropic thermopower of 
unstrained (strained) Bi2Te3 and Sb2Te3 at room tem- 
perature is shown, respectively. Blue thick (red thin) 
solid lines represent the in-plane thermopower 5*11 under 
electron doping (hole doping). The corresponding cross- 
plane thermopower S*^ is shown as a dashed line. The 
black dashed-dotted lines in fig. 5(d) emphasize the ex- 
pected doping dependent behaviour of the thermopower 
for parabolic bands, following the PlSARENKO-relation^^. 
For both tellurides we found, that the anisotropy of the 
thermopower shows a weak dependence on the strain 
state. However, for strained Bi2Te3 (see fig. 5(b)) the 
thermopower anisotropy under hole doping almost van- 
ishes, leading to S^^ ^ S^. It is worth noting, that the 
anisotropy of the thermopower is less pronounced for hole 
doping, than for electron doping for Bi2Te3 and Sb2Te3 in 
both strain states. As shown by the black dashed-dotted 
lines in fig. 5(d), the dependency of the thermopower 
on the charge carrier concentration differs from the Pis- 
ARENKO-relation^^ under sufficient high electron doping. 
This indicates, that the nonparabolicity of the energy 
bands has a noticeable impact in the investigated doping 
regime and should not be omitted by applying parabolic 
band models. 

Actually, changes for the absolute values of the ther- 
mopower can be found for both telluride systems under 
applied strain. In fig. 6 the relative change for the in- 
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plane component S^^ for both tellurides under in-plane 
strain is given. To compare the changes with the lat- 
tice constant, we relate the in-plane thermopower S^^ at 
the smaller lattice constant asbTe to the value at the 
larger lattice constant asiTe for both compounds. The 
doping was fixed to = 1 x 10^^ cm^"^ for Bi2Te3 and 
= 1 X 1020cm-3 for SbaTcg as done for fig. 2. Fig- 
ure 6(a) shows, that in the relevant temperature range 
between 350K and 450K the thermopower increases for 
Bi2Te3 under compressive strain for both p and n dop- 
ing by about 15-20%. For Sb2Te3 a decrease is ex- 
pected under tensile strain at electron doping and nearly 
no change under hole doping (see Figure 6(b)). With 
nearly all values above 1 for Bi2Te3 , as well as for 
Sb2Te3, it is obvious, that higher values of the ther- 
mopower require a smaller unit cell volume. One can 
expect, that the volume decrease causes a larger density 
of states and thus a shift of the chemical potential to- 
wards the corresponding band edge, connected with an 
increase of the thermopower S. However Park et al}'^ re- 
ported an unexpected increase of 16% for the in-plane 
thermopower S^^ of Sb2Te3 under p-doping (T=:300K, 
N = 1.32 X 10^^ cm"'^) if the material is strained into 
the Bi2Te3 structure. In the same doping and tempera- 
ture regime we find a slight decrease of 4% for S^^ . 

Comprising the statements on the electrical conduc- 
tivity and the thermopower, the related powerfactor for 
both tellurides in their bulk lattice and in the strained 
state are compared in fig. 7. It is well known, that 
optimizing the powerfactor aS^ of a thermoelectric al- 
ways involves a compromise on the electrical conduc- 
tivity a and the thermopower S^^ . Due to the inter- 
dependence of <J and S it is not advisable to optimize 
the powerfactor by optimizing its parts. In Figure 7(a) 
and (d) the doping dependent anisotropic powerfactor of 
unstrained Bi2Te3 and Sb2Te3 at room temperature is 
shown, respectively. Blue thick (red thin) solid lines rep- 
resent the in-plane powerfactor PFj under electron dop- 
ing (hole doping). The corresponding cross-plane power- 
factor PF^ is shown dashed line. Under p-doping 
both unstrained materials show a maximum powerfac- 
tor near carrier concentrations of ~ 4 x lO^^cm"'^. 
Absolute values of 35^W/cmK^ and 33/iW/cmK^ were 
found for unstrained Bi2Te3 and Sb2Te3, respectively, 
which is in good agreement to experimental and theoret- 
ical findings^'^"*'^^. Under electron doping the absolute 
values of PF|| (thick blue lines in fig. 7) were found to be 
distinctly smaller. This is due to smaller absolute val- 
ues of the thermopower for electron doping compared 
to hole doping (see Figure 2) and apparently smaller 
in-plane electrical conductivities a^^ at fixed carrier con- 
centrations. As a result, a powerfactor of 18yuW/cmK^ 
and 8 //W/cmK^ can be stated for unstrained Bi2Te3 and 
Sb2Tc3, respectively, under optimal electron doping. We 
notice, that the powerfactor for unstrained Sb2Tc3 is 
monotonically increasing for electron carrier concentra- 
tions of 6 . . . 30 X lO^^cm"^. This behaviour can 
be linked to a deviation of S*,, from the Pisarenko- 
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FIG. 7: (color online) In-plane (solid lines) and cross- 
plane (dashed lines) doping-dependent powerfactor at 300K 
for (a) Bi2Te3 in the Bi2Te3 structure, (b) Bi2Te3 in the 
Sb2Te3 structure, (c) Sb2Te3 in the Bi2Te3 structure and (d) 
Sb2Te3 in the Sb2Te3 structure. Electron (hole) doping is 
presented as blue thick (red thin) line. The charge carrier 
concentration is stated in units of e/uc (1/cm^) at the bot- 
tom (top) X-axis. 



relation under electron doping. While it is expected, that 
the thermopower will decrease for increasing carrier con- 
centration, S*!! was found to be almost constant in an 
electron doping range of A^ ~ 6 . . . 30 x 10^^ cm""^ (see 
fig. 5(d)). For the investigated electron doping range of 
A^ - 6 . . . 30 X 10^^ cm"3 the chemical potential ^ at 300K 
is located around 300 . . . 450 meV above the VBM. As can 
be seen from the band structure for unstrained Sb2Te3 in 
fig. 1(b) (black, solid lines) fiat non-parabolic bands near 
the high symmetry point Z dominate in this energy region 
and most likely lead to an increased thermopower. This 
feature is more pronounced for unstrained Sb2Te3, than 
for strained Sb2Te3 (red, dashed lines in fig. 1(b)). Sim- 
ilar statements can be done for strained and unstrained 
Bi2Te3 (see fig. 1(a)). We note, even though this picture 
is convincing, it is difficult to link such specific anoma- 
lies to the band structure on high symmetry lines, as the 
underlying TDF is an integral quantity over all occupied 
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states in the BZ. 

Under applied in-plane compressive strain for 
Bi2Te3 (ref. Figure 7(b)) and tensile strain for 
Sb2Te3 (ref. Figure 7(c)) the obtained changes in 
the powerfactor are noticeable different for both tel- 
lurides. While for Bi2Tc3 a decrease of the maximal 
powerfactor PF|| of about 27% and 23% for n-doping 
and p-doping was found, the strain shows nearly no 
influence on the powerfactor for Sb2Te3. At a carrier 
concentration of about ^ 3 x 10^^ cm~'^ the decrease 
in PF|| for Bi2Te3 is about 17% and 28% for n- and 
p-doping, respectively, while in the work of Park et 
al.^'^ a slight increase of PF|| under strain and hole 
doping is reported. Obviously this tendency has to be 
understood by analyzing the constituent parts and 
S'li . For comprcssivcly strained 612X03 at a hole carrier 
concentration of about iV ^ 3 x 10^^ cm^'^ the electrical 
conductivity decreases by about 39% to 330 (f2cms)~^. 
At the same time S^^ increases by about 9%, as shown 
in Figure 6(a). This results in the overall decrease 
of about 28% for PF|| . Under electron doping of 
A^ '-^ 3 X 10^^ cm~'^ no influence of strain could be found 
for S'li at room temperature (see solid blue lines in Figure 
6(a)). Thus, the decrease of PF|| under electron doping 
can be largely related to a decrease of the electrical 
conductivity under applied compressive strain. By 
detailed evaluation of the effective mass eigenvalues and 
eigenvectors we found a decrease of about 15% for the 
in-plane electrical conductivity of Bi2Te3 under applied 
strain in the low-temperature and low-doping limit'^^^'^^. 
The discussion can be made in the same manner for 
Sb2Te3^^''^^. The fact, that strain-induced effects in 
a and S tend to compensate each other was already 
reported for the case of silicon^^. 

As mentioned before (summarized in fig. 3 and fig. 4), 
we found a strong anisotropy in the electrical conduc- 
tivity with '^i\/a± 3> 1. The clearly preferred in-plane 
transport in both bulk tellurides is also reflected in the 
cross-plane powerfactor PF^^ (dashed lines in Figure 7), 
which is clearly suppressed for all strain states. It is ob- 
vious that PF^ is more suppressed for electron-, than for 
hole-doping. 

Nonetheless, we want to include experimental find- 
ings for the thermal conductivity to our calculations, to 
give an estimation for the figure of merit ZT in-plane 
and cross-plane. In Ref. 49 K|| = 2.2W/mK, k± = 
l.OW/mK, and K|| = 7.5W/mK, k± = 1.6W/mK for 
unstrained Bi2Te3 and Sb2Te3 are given, respectively. 
With this we find maximal values for the figure of merit at 
room temperature and optimal hole doping of ZT^^ ^ 0.48 
and ZTj_ - 0.41 for unstrained Bi2Tc3 and ZT\\ - 0.13 



and ZTj_ ^ 0.23 for unstrained Sb2Te3. We note, that 
the figure of merit ZT maximizes at slightly lower carrier 
concentration than the powerfactor aS'^ shown in fig. 7. 
This can be linked directly to an increasing electronic 
part of the thermal conductivity k^i with increasing car- 
rier concentration^^'^'^. 

V. CONCLUSION 

In the present paper the influence of in-plane strain 
on the thermoelectric transport properties of Bi2Te3 and 
Sb2Te3 is investigated. A focussed view on the influence 
of strain on the anisotropy of the electrical conductiv- 
ity (T, tlicrmopower S and the related powerfactor aS'^ 
could help to understand in-plane and cross-plane ther- 
moelectric transport in nanostructured Bi2Te3/Sb2Te3- 
superlattices. Based on detailed ab initio calculations 
wc focussed mainly on band structure effects and their 
influence on the thermoelectric transport. For both 
tellurides no reasonable decrease of the anisotropy for 
a and S could be found under strain, while in prin- 
ciple the anisotropy for a and S is more pronounced 
under electron doping, than at hole doping. Thus a 
favoured thermoelectric transport along the z-direction 
of Bi2Te3/Sb2Tc3-hcterostructures due to superlattice- 
induced in-plane strain effects can be ruled out and a 
clear preference of p-type thermoelectric transport can 
be stated for Bi2Te3 and Sb2Te3 and their related epitax- 
ial heterostructures. The absolute value of the in-plane 
thcrmopower S'li was increased under reduced cell vol- 
ume, which is in contrast to recent findings by Park et 
al}\ 

We found, that even if thcrmopower or electrical con- 
ductivity are enhanced or decreased via applied strain, 
they tend to compensate each other suppressing more 
distinct changes of the powerfactor under strain. We 
found the thermoclcctrically optimal doping to be in the 
range of A^ ~ 3 . . . 6 x 10^^ cm^'^ for all considered sys- 
tems. Our assumption of an anisotropic relaxation time 
for Bi2Te3 states that already in the single crystalline 
system strong anisotropic scattering effects should play 
a role. 
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